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Abstract. Gaussian and Chiral /^-Ensembles, which generalise well known orthogonal (/3 = 
1), unitary (/3 = 2), and symplectic (/3 = 4) ensembles of random Hermitian matrices, are 
considered. Averages are shown to satisfy duality relations like {/3,N,n} <4> {4//3,n, N} 
for all /3 > 0, where N and n respectively denote the number of eigenvalues and products 
of characteristic polynomials. At the edge of the spectrum, matrix integrals of the Airy 
(Kontsevich) type are obtained. Consequences on the integral representation of the multiple 
orthogonal polynomials and the partition function of the formal one-matrix model are also 
discussed. Proofs rely on the theory of multivariate symmetric polynomials, especially Jack 
polynomials. 
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1. Introduction 

It is the purpose of this article to obtain new dualities between different ensembles of 
random variables inspired by matrix models. Simple duality relations are given below, after 
a short review of the /3-Ensembles. 

1.1. /3-Ensembles. Let x = (x±, . . . , x?{) denote a set of A" random variables. Moreover, let 
(3 be a non-negative real number; it is called the Dyson index. The joint probability density 
function (p.d.f.) for the Gaussian /3-Ensemble (G/3E) is 

1 N 
P G /3E N (x)dx = | [ \xi -xjf ]Je~^dxi. (1.1) 

G/3Ejv l<i<j<iV i=l 

The average of a function / over the ensemble is given by 

(/(*)) ^ ■= I f(x)P G pE N (x)dx. (1.2) 

Note that Zq^ n is chosen such that the average of the identity equals one. Physically, the 
density Pg/3E n can be interpreted as the Boltzmann factor of a classical (two-dimensional) 
Coulomb gas at temperature 1/(3 |17j . There is also a quantum mechanical interpretation 
of the density: ^PqpEn is the ground state of a iV-body problem with pairwise potential 
interaction of the form 1/r 2 [9]. The chiral Gaussian /3-Ensemble (chflE) is defined similarly: 



P ch/3El(x)dx = — — Yl \xi-Xjff\x7e Xi dxi, (1.3) 

ctpE N l<i<j<AT i=l 

where the variables Xi are positive reals and the real part of 7 is greater then —1. Note that 
the Gaussian and chiral Gaussian ensembles are sometimes called Hermite and Laguerre (or 
even Wishart) ensembles respectively (cf. [11} I15j). 

In this article, the /^-Ensembles are considered from a Random Matrix Theory [21} [M] 
perspective. Let X = [Xij] be a random N x N Hermitian matrix whose entries are real 
(/3 = 1), complex ((3 = 2), or quaternion real ((3 = 4) (see Appendix A for more detail). The 
trivial (3 = case corresponds to a real diagonal matrix. It is a classical result that Eq. (|l.ip 
provides the joint density for the eigenvalues of a Hermitian matrix X when the latter is 
drawn with probability 

e" trX2 ((iX), (1.4) 

where (dX.) stands for the normalised product of all the real independent elements of [dX,^j\. 
Similarly, if a positive definite Hermitian matrix X is distributed according to 

(detX)^e- trX ((iX), (1.5) 

then the p.d.f. of its eigenvalues is given by Eq. (jl.3p . Moreover, set X = Y^Y, where Y is 
aJVi x N2 rectangular matrix such that N\ > N2. One can show that if Y has a Gaussian 
distribution, then the p.d.f. of X is given by JT3J) with 7 = (/3/2)(Aq - N 2 + 1 - 2/(3) [2H 
Chapter I]. 

It is worth mentioning that Eqs. (11. ip and (11. 3p can also be realised, for all (3 > 0, as the 
eigenvalue p.d.f. of tri-diagonal real symmetric matrices [15J. In the large A^ limit, scaled 
versions of these tri-diagonal matrices can be seen as stochastic differential operators 
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1.2. Some results. The aim of this paper is to prove the equivalence of some averages over 
/3-ensembles and other averages over 4//3-ensembles The averages contain a source matrix 
S of size n, and an external field matrix F of size N. These matrices are Hermitian with 
eigenvalues s = (si, . . . ,sn) and / = (aci, . . . , rcjv) respectively. Roughly speaking, the duality 
relations obtained here read 



f 1 




■ A/ (3 ' 


TV 




n 


< F 




S > 


n 




AT 


S 

V J 




F 



In order to formulate the results precisely, further notation has to be introduced. Let 1 be 
the identity matrix and let 

(FWxeGPE = J (^X)e- trX2 F(X)/ J (dX)e~ trX2 , (1.7) 

which is the average of the complex-valued function F over the Gaussian /3-Ensemble of 
matrices with p.d.f. (|1.4p (not to be confused with the average over the eigenvalues). Special 
realisations of Eq. fll.6p are given in the following propositions. Only expression relative to 
the Gaussian ensembles will be displayed for the moment; dualities for the Chiral ensembles 
will be given in Sections 3 and 4. 

Proposition 1. Suppose that (3 = 1,2 or 4 and that (3' = 4/(3. Then 

e-Wndet^lii^xV"-) 

e~^ 2 (f[ det (Y ± iy^/.l) e 2trYS \ . (1.8) 

The previous result concerns the average of products of characteristic polynomials. Special 
cases of this duality previously appeared in the literature. For instance, the duality for 
moments of characteristic polynomials, which corresponds to s% = . . . = s n and F = 0, has 
been proved for (3 = 1, 2, 4 and conjectured for all j3 by Mehta and Normand For F = 
(i.e., no external field) and (3 = 2, but for distinct s, this has been observed by Fyodorov and 
Strahov in [23]. Based on a previous work [6] and the supersymmetric method, Brezin and 
Hikami [7] have obtained dualities reproducing Eq. (|1.8j) in the case F = and (3 = 1,2,4. 
The problem with a non-zero external field and (3 = 2, has been recently solved by same 
authors [8]. 

A less common duality of the type (|1.6p . but not affecting (3, concerns products of inverse 
characteristic polynomials. 
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Proposition 2. Suppose that /3 = 1,2 or 4 and that (3' = 4//3. Assume moreover that the 
variables s and f have non-zero imaginary parts. Then 



More general dualities will be exposed in Sections 3 and 4 after the introduction of multi- 
variate hypergeometric functions in Section 2. To the best of the author's knowledge, explicit 
dualities valid for (5 general and N finite, first appeared in a 1997 paper by Baker and For- 
rester [3]. For instance, Eq. 5.31 in [3] is equivalent to Eq. (jl.8p above with F = and 
s\ = . . . = s n = t say, and (3 integer. The latter expression has been used in [11] for calculat- 
ing asymptotic corrections to the global eigenvalue density when (3 is even. The present work 
can be considered, to some extend, as a continuation of |3j. Note also that Proposition 7, 
for s\ = . . . = s n = t and / = 0, provides a proof of the formula conjectured by Mehta and 
Normand for all (3 > [33 Eq. 3.29]. 

Closely related are the "particle-hole" dualities observed, when (3 is rational, in the limit 
N — > oo of the correlation functions for the Sutherland model or equivalently, for the Circular 
/3-Ensembles (see for instance [39J and references therein). This can also be interpreted as a 
strong- weak coupling duality (see for example [lj). 

Inter-relationships between orthogonal, unitary and symplectic ensembles have been also 
considered, at the level of the joint eigenvalue p.d.f. itself, for eigenvalues respecting in- 
terlacing inequalities [22J. Very recently, this work has been generalised by Forrester in 
[23]. The duality found in the latter reference can be summarised as follows : setting 
(3 1 = 2(r + 1) = 4//3, the joint distribution of every (r + l)-st eigenvalue in a /3-ensemble 
is equal to that of a /^'-ensemble if the eigenvalues are properly ordered. 

Before going further into the study of dualities, a last comment is in order. For (3 = 
1,2 and 4, expectation values of products and ratios of characteristic polynomials can be 
expressed in terms of determinants or Pfaffians (see the extensive study by Borodin and 
Strahov [5] and references therein). The size of the determinants depends on the number n of 
characteristic polynomials but not on the size N of the random matrix. Thus, one can exploit 
the determinant formulae for calculating the asymptotic behaviour of the correlation functions 
when N — * oo. For general /3, no such determinantal formulae are available. Then, solving 
a /3-matrix model amounts to finding a reduced integral representation for the correlation 
functions. In other words, the aim is to obtain integral formulae whose dimension does not 
depend on the size of the matrix, thus providing a representation for the correlation functions 
that, in theory, allows to take the limit N — > oo. 





'n 



(1.9) 



2. Preliminary definitions 



This section furnishes a brief introduction to the theory of symmetric polynomials. More 
detail can be found in [2H l26l \32[ [43] 
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2.1. Partitions. Let A = (Ai, A2, . . . , ) denote a partition of n; that is, a sequence of non- 
negative integers such that 

Ai>A 2 >...>0, |A|:=^Ai = n. (2.1) 

i 

One usually writes (Ai, . . . , A^, 0, . . . , 0) = (Ai, . . . , A^), where t = £(X) gives the number of 
non-zero parts in A. A partition can also be expressed as follows: A = (. . . , 3™ 3 , 2 n2 , l ni ) 
where = nj,(A) is the number of parts of A that are equal to k. The conjugate partition of 
A, written A' = (A^, X' 2 , ■ ■ .), is such that 

A^ = #{Ai e A : Aj > A;} (2.2) 

To each partition, we associate a diagram by drawing Ai boxes on the first row, then A2 
boxes under the first row, and so on, all boxes left justified. See Fig. [TJ The conjugation 
of the partition then corresponds to the transposition (as for matrices) of the diagram. The 
arm of a point s = in A, written a\(s), is the number of boxes to the right of in 
the ith row of the diagram A. As illustrated in Fig. O similar definitions exist for the leg I, 
co-arm a', and co-leg /' of the point s = 



Figure 1. Diagram and conjugate diagram 



A' 



The dominance ordering of partitions is defined as follows: 

k 

\>/j <^ ^(Ai - m) > VA;. (2.3) 
i=i 

This ordering is partial. Note the obvious property: 

A > fi <*=> A' < //. (2.4) 



Four functions on partitions will be frequently used in the article. They are: the generalised 
factorial 

,_- rT r(n-(j-l)/q + A i ) - 

-ii Y(u-{j-l)/a) ' (2 ' 5) 

the specialisation coefficient 



b[ a ' N) = l[(N-l> x (s) + aa> x (s)); (2.6) 
sex 
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and the "lower" and "upper" hook-lengths [43J Q 

Y[ (h(s) + 1 + aa\(s)) and = (l\(s) + a + a a\(s)) ; 



(2.7) 



sex 



sex 



Note that b[ a ' N) =a\ x \[N/a]^ N) . 



Figure 2. Arms and legs of s = G A 



l'x(s) 



lx(s) 



2.2. Jack polynomials. Let x stand for the ordered set (x\,...,xn). A function f is 
symmetric if it is invariant under permutations of the variables; that is, f(x) = f(x a ) for 
any element a of the symmetric group Sn, where x a = (x^m, • • • , %a(N))- The set of all 
symmetric polynomials in iV variables, whose coefficients are rational functions of a, form 
an algebra over the ring Q(a), witten P Sn . A symmetric polynomial that is homogeneous of 
degree n can be decomposed into the monomial basis {m^ : |A| = n}, where 

and n(A)! := m(A)!n 2 (A)! • • • . 

Another important basis for P Sn is provided by products of power sums {p\}, where 

N 

Px{x) ■.= px 1 (x)---px l (x), Pk ( x ) = J2 x i- ( 2 - 9 ) 

i=i 

The combinatorial (or Fock space) scalar product in the algebra of symmetric polynomials 
can be defined by 

((px\p,)) ia] =a^z x 5 x ^ (2.10) 
where z\ = \\ k k nk k\. The parameter a is related to Dyson's (3 index as follows: 

2 

(2.11) 







The Macdonald automorphism will play an important role in the following paragraphs; it 
is given by 

UJkPn = {-iT^kpn (2.12) 



h} a) and h " are respectively noted h\ and d' x in |21] , while they are noted C\ and c' x in 
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and satisfies 

(M\g)) {a) = ((^ /a f\g)) ia) = «<^/b» (1) - (2-13) 

for any symmetric polynomials / and g. 

The (monic) Jack polynomials, denoted by P\(x) = P\(x;a) or by P^*'(x), generalise 
many important symmetric polynomials: 



P (») 
A 



r e\>, 


a = 





s\, 


a = 


1 




a = 


2 




a = 


oc 



(2.14) 



where e\, s\, and Z\ respectively stand for the elementary, Schur, and Zonal polynomials. 
A 



-fff^ is the unique symmetric polynomial with coefficients in Q(a) satisfying 



(1) P x = mx + y^ c\^,(a)mx (triangularity) 

n<\ (2.15) 

(2) ((Px\P^)) {a) = ||^ a) || 2 V (orthogonality) 

Alternatively, the Jack polynomials can be considered as the unique triangular symmetric 
polynomials that comply with 

D P\(x) = e\(a)P\(x) (eigenfunctions), (2-16) 

where 

^ = E(-?C+2£^A (2.17) 

i=l \ j# X * X * j 

The eigenvalues are given by e\(a) = Y2i>i ( a (^ ~~ ~ (* ~~ The following speciali- 

sation and normalisation formulae will be used later: 

(a) (a) 

The Jack polynomials possess several remarkable properties. Amongst them, it is worth 
mentioning their duality (see eq. (|2.13p ) 

UaP{ a) = iff with «':=-, (2.19) 

A ||p( a )||2 A a 



tttt — —— = — i_pj a )( z ) p A (a) ( y ). (2.20) 

1111 (1 -X,-Vi) 1 / Q ^||p( a )||2 V ^ 



and the Cauchy type formula 

AT M 

i=lj=l ^ ly J> A I^A II 

By applying the former property to the latter equation, one gets 

N M 

]J II( 1 + = E p v (2-21) 

i=lj=l A 
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2.3. Multivariate hypergeometric functions. In multivariate analysis, the Jack poly- 
nomial P\(xi, X2, ■ ■ •) plays a role that is similar to that played by the monomial for 
functions in one variable. As an example, the multivariate hypergeometric functions of two 
sets of variables are given by [211 144] : 

p F q ( a '(ai,.. . ,Op]bi,.. . ,b q ;xi, . . . ,x N ;yt, . . . ,yjv) = 

V « |A| [aii a) ---K][ a) P { x a) (-)P ( x a) (y) (2 22) 

There is no simple "explicit forms" for the latter functions in general. However, one can 
prove that [27] 

JV 

T ^( Xl ,...,x N ;t,...,t) =l[e tx > (2.23) 

i=l 

N 

iT ^(a;xi, . . . ,x N ;t, . . . ,t) = Y[{1 - txi)- a . (2.24) 

8=1 

Softwares are also available for computing "truncated" hypergeometric functions; they are 
based on the exact calculation of Jack polynomials (e.g., see [28 J. 

As explained in [3], the p J-q^ a ^s provide generating functions for the multivariate Hermite 
and Laguerre polynomials, respectively written and C x a ^\ These polynomials were first 
defined by Lassalle in [30], [31] . Explicitly, 

e -P.(y) ^ W(2x;y) = £ J -l^ n[ a \x)P x a \y), (2.25) 



where 



M»,N) = ^h^Pi a \l N ), (2.26) 
e -Mv) ^H (7 + q . x . y) = Y j —J^c[ a ^\x)Pi a \y), (2.27) 

B x (a,j,N) = 2^[ 1 + q}[ a) A x (a,N), q = 1 + (JV - l)/a. (2.28) 

Multivariate orthogonal polynomials are deeply connected to the Dunkl operators and Calogero 
systems [3] Q2] [38]. From these relations, it is possible to show, for instance, that the Her- 
mite and Laguerre polynomials in many variables provide orthogonal bases for the algebra of 
symmetric polynomials: 



and 



where 



(h[ 2/ V(x)H^(x)) = A x {2/(5,N)5 x ,» (2.29) 

(c[ 2/ ^\x)C^(x)) = B X (2/^,N)5 X ,,. (2.30) 



2 A11 polynomials involved in Eqs. (I2.25|) and (I2.27P are monic while those in [3] are not. One goes from 
a convention to another via H x a) (x) = 2^ x ^ P { x a) (l N )H x (x; a) and ^"'^(x) = (-1)I A I \\\\P { x a) (l N )Ll(x; a). 
Additionally, P { x a) (x) = a~^h x a) C{ a) (x) /\\\!. 



DUALITY IN RANDOM MATRIX ENSEMBLES 



9 



3. Expectation values of Jack polynomials 

In this section, it is proved that the averages of Jack polynomials over the /3-ensembles 
defined in (jl.ip and fjl .3 j) enjoy simple duality properties. Recall that the Jack polynomials 
form a basis for symmetric polynomials. Thus, in principle at least, dualities involving Jack 
polynomials provide general tools for reducing the dimension of the correlation functions. 

3.1. Gaussian ensembles. The starting point is a conjecture by Goulden and Jackson [25] 
first proved by Okounkov in [36]. Another proof is given below. 

Lemma 3. Let X be a partition of an even integer and let fi be the partition (2l A l/ 2 ). Also, 
let {{f\g)) denote the combinatorial scalar product denied in Eq. (|2.1Qj) . Then 



Alternatively, 



pW\x)) = b^ N) coeSP^\ (3.2) 

'xGQSEjv Pm 



where coeff/ denote the coefficient of x in f . 

X 

Proof. It is simpler to set a = 2/j3. The orthogonality (|2.29[) and Eq. (|2.25p immediately 
imply [3, Proposition 3.8] 

( ^o (a) (2x;yW a) (2*;z)) = e^ + ^ z \^ a \2y; z). (3.3) 

\ / x£G/3E 

Then, setting z = in the last equation and recalling that -Pa(O) is equal to 1 if A = (0) and 
otherwise, one gets 

( ^o (a) (2x;y)\ = e^ . (3.4) 

Now, the left-hand-side is developted in terms of the Jack polynomials and the combinatorial 
scalar product with respect to P^ is taken (see Eq. (|2.15|) ). This leads to 

|A|/2 



(a) A V / 

which, by virtue of Eq. (|2.18p . is equivalent to the first equation of the lemma. The second 
equation simply follows from 

«^|p M » (Q) = (2a)l A l/ 2 (|A|/2)! if ^ = (2' A I/ 2 ). (3.6) 

□ 

Proposition 4. Set x' = ix^/2/P and /3' = 4//3 > 0. Then, for every N > £(X) and 
N' > i(X'), 

I \ I ^!V) \ f37) 

VA v II x6G/3E]v \ r \' \ L II x£G/3'E N , 
Note that, as a consequence of the orthogonality (|2.29j) . the expectation values are non-zero 
only if \X\ is even. 
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Proof. Assume N > £(X) and |A| even. Following the notation used for the previous lemma, 
one has 

The Macdonald automorphism (|2.12p is then exploited for getting 
Moreover, the use of Eq. f)2. 13|) yields 

= (- 1 ) <M2<M ^|« w » J 'i° ) W (1/< "'- (3-10) 

But according to the duality Eq. (|2.19p of the Jack polynomials, the last equation can be 
rewritten as 

(^SL, E „ = ( - i) ' M2 ' < " wM - <3ji) 

where a' = 1/a. Note that the last displayed equation is non-zero only if N' > where 
iV 7 denotes the number of variables in pj? \ Besides, from the definition (|2.7p of the lower 
and upper hook-lengths, one easily obtains 

h X a)=a \X\ h («') =a 2lM h (-') (3 . 12) 



Hence, 

(iSfk) =(-l)'""2'«a 2 '«^«<»W) < "' ) . (3.13) 

Finally, Lemma 1 or Eq. (|3.8|) is used once again to obtain 

\- r X ^ n xeG0E N » A' ^ ^ / yGG/3'Ejv, 

which is the desired result. □ 

The simplest example of such a duality is certainly the average of moments of the charac- 
teristic polynomial (cf. [35]). Indeed, from 

(detX) n = (si • ..x N ) n = P { ffi{xx, ...,x N ) (3.15) 

and the substitution of the the latter formula into Eq. (|3.7p . one concludes that 

<(detXr> XeG/3Ejv = (-2/(3)^ <(detX)^) XeG/3 , En , (3.16) 

for = 4//3' = 1, 2, 4. More generally, for all j3 > 0, 

<(*i • • • ^T) xeGmN = (-2//?)^ 2 ((X! • • • x n ) N ) x ^, En . (3.17) 
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Another simple manifestation of Proposition is related to the Hermite polynomials. It is 
well know that an orthogonal polynomial of degree N has a N x N matrix integral represen- 
tation. However, a classical orthogonal polynomial can also be realised as a single integral. 
For instance, the Hermite polynomial of degree N can be written either as 

N 

H N (t) = 2 N (T[(t-x l )) =2 A Ydet(tl-X)\ (3.18) 

X 1 - 1 / XGG2E \ /XGG2E 

i=l 



or as 



2iv /-00 

H N (t) = —= dxe~ x (t±ix) 

V n J-oo 



N 



(3.19) 



Proposition [4] allows to go directly from one representation to the other without referring to 
the orthogonality, determinantal representation, or differential equation. On the one hand, 



N N N 



Y[(t -x i )=Y J t N - n (-l) n e n ( Xl , ...,x N )=J2 t N - n (-irP$ } ( Xl , x N ), (3.20) 

i=l n=0 n=0 



where a = 2/(3, and 



pfrki N ) = r ), (3-2i) 



so that 

2=1 ™=° V 7 " G pE N 

On the other hand, Proposition 1 and (l n )' = (n) immediately imply that the last equation 
is equal to 



n— u \ ( n ) V / / Qg; El n— u 

(3.23) 



Consequently, for all /3 = 2/a > 

N roc / /— \ N 



IIC*-^)) «™ dxe- x \t-\^x) N =(^) H N {t/^) (3.24) 

3.2. Chiral ensembles. Lemma [3] and Proposition 0] can be easily adapted to the chiral 



case. 



Lemma 5. Let fi stand for the partition (l' A ') ; i.e., £(n) = |A|. Then, for all partition X, 
(3 > 0, and 7 > — 1, 



where q = 1 + /3(JV — l)/2. Equivalently, 

P^\x)) =[ 7 + g] f«6?/^ooeffP^. (3.26) 
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Proof. The orthogonality (j2.30j) and Eq. (12.270 give [31 Proposition 4.11] 

( o^i (Q) (7 + q; x; y) o^i (Q) (7 + 5; x; z)) = ^ ( 7 + q- y- z) (3.27) 

for a = 2/ p. Then, one sets z = and proceeds as in Lemma El □ 

Employing the same method as that exposed in the proof of Proposition [H one shows 
that the last lemma and the duality (12.190 for the Jack polynomials imply the following 
proposition. 

Proposition 6. Set 0' = 4//9 > 0, q = 1 + (JV - l)/ty2, and q' = 1 + (N' - 1)2/(3. Moreover, 
assume that N > i(X), N' > £(X'), 7 > -I, and i > -I. Then, 



\pWP)(lN\ U/ + /l( 2 //3') \p(Vl3') (1 N')/ 

As an example, choose A = (n ) and 7' = 2(7 + l)//3 — 1. Direct manipulations then give 

h + Ql x I P 



/l(2//3') 



(3.29) 



This implies that moments of the determinant also satisfy a simple duality in the Chiral 
/3-Ensemble: 

<(*i---^rW; = (£f (^■■■*n) N ) xecW< • (3-30) 

It is worth mentioning that the average in the Gaussian ensembles are limit cases of similar 
averages in the Chiral ensembles. This can be understood as follows. Set y = 7 + \/2rjx with 
7 > 0. Trivial manipulations give 

7 - 1 eVe- ? ' = e~ x2 + 0( 7 ~ 1/2 ) (3.31) 
uniformly when 7 goes to 00. Then, simple changes of variable allows one to conclude that 

Km (F(-V^ + y/V^)) yech0El = <*X*) W w (3-32) 

for all multivariate polynomials F, positive integers iV and (5 > 0. For instance, by setting 
7' = 7 and taking the above limit in both sides of Proposition (U one readily establishes 
Proposition HI 



4. Source - external field dualities 



In this section, averages of (inverse) characteristic polynomials in /3-Ensembles with exter- 
nal fields are studied. Central to the analysis is the use of Dunkl transforms (see for instance 
PI). 
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4.1. Gaussian ensembles. The following formula is a generalisation of the Fourier trans- 
formation O [38] : 

e~ P2(y) / ^"o (2//3) (2x;y)F(x)\ = e^'^ F{y) (4.1) 

\ / xeGfSE 

where 

A (q) _ ^ ( d | 2 ^ 1 ^ A. (4.2) 

I 9a;,- 

i=l \ j+i J J 

The equation holds true for any analytic function F. It can be seen as a consequence of Eq. 
(|3.2p and the Lassalle formula 

«W(a;) = e -^ o) pW(i). (4.3) 

Now imagine that there exists a function G[s' 7 /*) of two sets of variables, 5 — (^l? • • • ; ^n) 
and / = (/i, . . . , /tv), satisfying for some a' 

A( a) G( S ;/)=A(«')G( S ;/). (4.4) 
Due to the commutativity of Ay and A^" \ this means that 

l.W 1 . (a) 

e* A f G(s-J)=e* A ° G(s;f). (4.5) 
Therefore, as a direct consequence of the latter formula and Eq. (|4,ip . 

e-»W ( ^o (2//3) (2x; x)\ = e -»W ( o^o (2//?,) (2y; /)G(y; /)) (4.6) 
where it is understood that j3' = 2/a' . 

Proposition 7. Set (3' = 4//3. Then, for all positive integers n and N, and all a > 0, 



g-P2(/) 



(4.7) 



yeG/3'E„ 

J/ m addition the variables s and f are not real, then 

I n N \ 

e -Mf) / JJ JJ fciafc)-^ o^o (2//3) (x,2/)\ 

I n N \ 

^ w nn^^r^^^^) . (4.8) 

Remark that in the last equation, (3 is not affected by the duality transformation. 
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Proof. Recall that a = 2/(3 and let 

n N 
3=1 k=l 

Following the above discussion on the Dunkl transform, it is sufficient prove that 

A^H(s;af) b = A^H(s;af) b (4.10) 
for appropriate values of a, b, and a' . On the one hand, by using 

4n^/)'=-«'E^n^/) 1 («d 

and 

?S«i-a/fc/fe-/i = 2^^, *i - afr Si - aft (4 ' 12) 

% k+l i k+l 

it is simple to show that 

On the other hand, exploiting 

^ u s i ~ a fk Si- sj 2 ^ ^ Sj - a/fc - o/fc ' i^ 15 ) 

i^tj A; J j^j k J 



and 



one finds 



Tr* Si - afk V s « ~ 

2,K 



'S^s + '£^) n( ^ ,4 - 16) 



Note that s and / are kept generic. Thus, imposing the equality of Eqs. (|4.13p and (|4.16p 
requires 

a 2 (b- 1) = (b- 1), a'a 2 6 = -l, a 2 = -ba. (4.17) 
The only solutions to the latter system of equations are either 

a = ±iy/a, 6 = 1, a' = 1/a. (4.18) 

or 

o = ±l, 6 = -1/a, a' = a. (4.19) 
This completes the proof of the proposition. □ 



DUALITY IN RANDOM MATRIX ENSEMBLES 15 

Consider for example the expectation value of a product of characteristic polynomials 
without external field, i.e., for / = 0. According to Proposition 

(nilf^V^)) =e P2{s) (h(iyj) N o^'H2y,-is)) . (4.20) 

\i =1 fc=1 ' xeG0E N \j =1 ' yeG/3'En 

The right-hand side can be written as 

e»« (pl^(iy) ^\2y,is)) yeG ^ . (4.21) 
Due to the Dunkl transform (|4.ip . the latter equation is also equal to exp ( — iAg ) P r ™?(s). 



4" s y x (N n ) 

Returning to the Lassalle formula ()4.3I) . it becomes clear that averages of products of char- 
acteristic polynomials are multivariate Hermite polynomials: 

<»l( s )=(nn( s i ± A/f xfc )/ • (4 - 22) 

This equivalence was first pointed out in [3] by considering a limit of the multivariate Jacobi 
polynomial. Another proof consists in showing that the expectation value satisfies a Calogero- 
like eigenvalue problem. 

Propositions Q] and [2] given in the Introduction are in fact special cases of Proposition \7\ 
The connection between the averages involving multivariate hypergeometric functions and 
the matrix integrals can be understood as follows. Let (3 = 1,2, or 4, and let (U'cHJ) stand 
for the normalised Haar measure for unitary matrices with real, complex or quaternion real 
elements respectively. It is a standard result [211 El] that 

1 N 

(dX) F(X) = - Yl dx Yl \ x i -Xjf(UU\J)F(UxU^). (4.23) 

i=l l<i<j<N 

for any Hermitian /3-matrix X = XJxTJ\ where x = diag(xi, . . . , scjv), and some constant C. 
When F is not invariant, i.e., when FiXJxXJ^) ^ F(x), the calculation of the expectation 
values containing F requires the use of the theory of zonal polynomials (see Chapter VII in 
Macdonald's classical book [32]). One can show for instance that, if X and Y are N x N 
Hermitian /^-matrices, 

(UtdU)pf //3) (XUYUt) = (4.24) 

In the last equation, p| 2 ^(XUYU^) = P^ 2 \z\, . . . , zn), where (z±, . . . , zn) denotes the 
eigenvalues of XUYU*. Now, by exploiting formula (|2.23p . one readily gets 



/' 



(UtdU)e tr ( XUYUt ) = ^o (2//3) (x;y), (4.25) 
and consequently, for any invariant function <?(X) = g(x) and = 1, 2 or 4, 

0.P0 (2//3) (*;/)<?(*)) nm} = 

/ (dX)e-^ 2 ^g(X)/J (dX)e" teX2 = <e trXF 5 (X)> XeG/3Ejv . (4.26) 
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Finally, the comparison of the latter formula with Eqs. (14.71) and (|4.8j) respectively establishes 
Propositions CD and EJ 

4.2. Chiral ensembles. There exists a Lassalle formula for the multivariate Laguerre poly- 
nomials. For a = 2/(3 > and 7 > — 1, it reads [3] 

C^\x)=e- A ^ ) P^(x), (4.27) 

where 

From formula (|4.2T|> and Eq. (13.27D , one can prove the following generalisation of the Hankel 
transform which is due to Dunkl [3j [38] : 

( ^i {2//3) (7 + q;x;-y)F(-x)) = e'^'^ F{y) (4.29) 

with and q = 1 + /3(iV — l)/2. The strategy for getting dualities in the Chiral /3-Ensemble is 
the same as in the Gaussian /3-Ensemble; that is, find a function G(s; f) such that 

A^ 7) G( S ;/) = A^')G( S ;/), (4.30) 

and exploit the Dunkl transform (|4.29j) in order to conclude that 

e^)( ^i (2//3) (7 + q;x;-f)G(s;-x)) = 

^1(S) /^/P')^ + q >.y._ s)G{ _y. f) \ (4 . 31) 

\ / yech/3'ET 

for some 0', 7', and q' = 1 + /3'(n — l)/2. 

Proposition 8. Set /3' = 4//3, 7' = 2(7 + — 1. Then, for all positive integers n and N , 
all a > and 7 > —1, 

n at 

a pi(f) 



UUU-h) 0^ 2m (7 + q;x,-f)\ 

(-l)^W (fl ft U - |/*) o-^V + </;</,"*)) • (4.32) 

Suppose moreover that s and f have non-zero imaginary parts, then 
I n N \ 

ePi(/) ( n n & + ^) _/3/2 o^i {2//3) (7 + ?; -/) ) 



V7 =1 fc=1 ' xech/JE^ 



I n N \ 

(-l)^V^) / J] n (yj + / fc )^/ 2 o^i (2//3) (7' + q';y, s) ) . (4.33) 

77ie latter equation is valid only if (3 > 2~f, which ensures 7' > —1. 
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Proof. Define 

n N 

ll(s;af) b = ilil(s3-af k ) b . (4.34) 

3=1 k=l 

Direct calculations lead to the conclusion that 

A /°' 7) W s a ff = W s a f? ( 4 - 35 ) 

for all n and N, if and only if, 

aba' = -1, a(b- 1) = (6- 1), a( 7 + 1) = -(7' + 6), a = -6a. (4.36) 

The proposition follows from the fact that the unique solutions to these equations are 

a = —a, 6=1, a = l/a, 7' = 0(7 + 1) — 1 (4.37) 

and 

a = l, 6 = -l/a, a' = a, Y = l/a - 7 - 1. (4.38) 

□ 

The comparison of Eqs. (I4.27p . (|4.29p . and ()4.32p leads to a simple generalisation of Eq. 
(I4.22p . that is, a random matrix representation of the multivariate Laguerre polynomial: 

^? ,) ( s )=(nn(^-^)) > 

\ j=l k=l V 7 / ^.gcl^^ 

where /3' = 4//3 and 7' = 2(7 + l)/3 - 1. 

5. Further applications 

5.1. Multiple polynomials. Ensembles of random Hermitian matrices with external field 
naturally lead to multiple polynomials when (3 = 2. These polynomials satisfy orthogonality 
conditions for more than one scalar products; they can be interpreted as multi-parameter de- 
formations of the usual orthogonal polynomials in one variable (see the review by van Assche 
in [2]). 

It is known that the multiple Hermite polynomials of type II can be defined as follows: 

H n (z) = e"3 trF2 (det(*l - X)e trXF > XeGUEjv (5.1) 

Here the polynomial is chosen monic of degree N while n = (nj., . . . ,nr>) is the multiplicity 
vector for the eigenvalues of F. In other words, F has rii eigenvalues equal to gi say, and 
Y2i n i = N. One can write for instance 

f = 9 n = {g?T--,9 n D D )- (5-2) 

As a consequence of Proposition 1, it is clear that the multiple Hermite polynomial also has 
the following realisation: 

( _-\N r J) 

H a (z) = ^—L- / dye^-^ 2 H(y - (5.3) 



3 The factor exp(-trF 2 /4) comes from the absence of F in the definition (|1.4|l of p.d.f. for GUE. See also 
Eq. (|A~9|l . 
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This formula has been first obtained by Bleher and Kuijlaars in [4] from the orthogonality 
relations satisfied by P n . 

There is another interesting application of Proposition 1. Indeed, the comparison of Eqs. 
(jl.8p and (|5.3p indicates that the multiple Hermite polynomials can be represented as (3- 
matrix integrals 



for (3 = 1, 2, 4, or as 



H n (z) = e-^(l[(z-^ Xj \ ^ Uf-f)) (5.5) 



2GG/3EJV 



\j=l 

for all /3 > 0. 

The multiple Hermite function of type I can also be written as a single integral [4]. 

Q n (z) = lim -9 / dye-(y- z+i ^ T\(y - fj/2) n '. (5.6) 

Note that taking the imaginary part is equivalent to closing the contour around the poles. 
Now suppose that (3/2 is a positive integer and that n« = (3rrn/2 for all i = 1, . . . ,D. In 
addition, construct a reduced sequence of variables (only the multiplicities are changed) 

f = g m = (gT 1 ,...,9D D )- (5.7) 

and set M = Y2j m j = (3. Then Proposition 2 furnishes new integral representations for 
the multiple Hermite functions of type II: 



Q n (z) = e -X/) Km -%lTUz-k-Xi)-V 2 ^W®(x'J)\ . (5.8) 

>0+ 7T \ ' 



M 



Similar representations can be deduced for the multiple Laguerre polynomials of type I 
and II by exploiting the results given in Section 4.2. 



5.2. Formal one-matrix models. One point about matrix models encountered in theoret- 
ical physics and combinatorics deserves to be clarified. 

Matrix integrals are often used for solving combinatorial problems related to 2D Quantum 
Gravity, such as counting graphs drawn on surfaces of a given genus (for a simple but precise 
introduction to the subject, see [45] : more advanced topics and recent references can be found 
in [H]). Questions of convergence of the integrals are not an issue when enumerating objects; 
only formal power series are considered. One typically looks at functions involving traces of a 
random matrix averaged over the Gaussian ensemble of Hermitian matrices (usually (3 = 2). 
Averages are obtained from the partition function of the one-matrix model: 
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For (3 = 2, Chekov and Makeenko [TO] have shown, by exploiting Schwinger-Dyson equa- 
tions, the equivalence of the one-matrix model to the Gaussian model with an external field. 
In fact, this remains true for (3 = 1 and 4. It can be proved quickly. First, set 

t x = -2u + i>trS _1 , t 2 = 2 + u 2 trS" 2 , and t k = v k tTS~ k for k > 3, (5.10) 



where u is a formal parameter, v = ±i v / 2//3, and S is a n x n Hermitian or anti-Hermitian 
whose eigenvalues s are non-zero (so, is essentially a kth power sum). Then, the formal 
development of ^ fc>1 t k tr X. k /k in terms of the eigenvalues of S yields 

expj-J^trXM = (detS)" iV ]Jdet(s i l -fX)exp (-trX 2 + 2utrX) . (5.11) 

Finally, from the substitution of the last equation in Eq. (|5.9p and the comparison with 
Proposition 1, it holds that the formal one-matrix model is equivalent to a Gaussian model 
with an external field. 

Proposition 9. Let X be a N x N Hermitian (3-matrix. Moreover, suppose that (3' = 4//? 
and that S is a n x n Hermitian 4/ ' (3-matrix satisfying Eq. (|5.10p . Then 

-trS 2 +JV« 2 / / [g \ N \ 

Z N ,p(h,t 2 ,...) = z N} p detgjV /det Y-i^l e 2trYS \ , (5.12) 

<U \ V / / YeG/3'E„ 

where 

z N3 = J (dX)e" trX2 = 7 r Ar / 2 (7r/2)^( JV - 1 )/ 4 . (5.13) 

The latter result remains true for all (3 > if one appropriately uses the generalised 
hypergeometric functions o-^o(2y;s) in place of e 2trYS (cf. Section 4.1). 

5.3. Matrices at the edge of the spectrum. Consider a random matrix belonging to 
the G/3Ej\r where (3 = 1,2 or 4. Suppose moreover that the size N of the matrix goes to 
infinity. Then, a classical result [20] says that that the eigenvalue correlation functions, when 
recentered and rescaled at the spectrum edge (i.e., Xi ~ V2N), can be expressed in terms of 
determinants (or Pfaffians) involving the Airy function and its derivative. Here it is proved 
that the average of products of characteristic polynomials can be expressed, at the edge of 
the spectrum, as a matrix Airy integral. For (3 = 2, this integral has been first introduced 
by Kontsevich in [29J when studying the asymptotic behaviour of the one-matrix model's 
partition function (15. 9|) . 



Let W and F be N x N Hermitian /^-matrices. The matrix integral of the Kontsevich type 
is defined as follows: 

Ai( 2 // 3 )(F) = a% /l3) J (dW)exp ^trW 3 + itrWF^) , (5.14) 

where aff = (27r)~~ N ~ N ( N ~ 1 ")/ a . Obviously, the standard Airy function is recovered for 
N = 1. Recall that the latter function satisfies Ai"(x) = xAi(x). From the derivation of 
Eq. (15.14p and the use of integration by parts, a simple generalisation of the Airy differential 
equation is obtained: 

A F Ai(F) = trFAi(F). (5.15) 
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In the last equation, Af stands for the Laplacian of the matrix F; that is, 

i=l Z a=0 1<i<j<N \°*vj 

In fact, the matrix integral (|5.14p is a special realisation of a more general multivariate 
Airy function (see Eq. (|4.25p ). 

Ai^(f)=aP I e^ w l ^ a \w;if) J] \w 3 -w k \ 2 ' a dw, (5.17) 

•^ RJV l<j<k<N 

which is valid for any a = 2/(3 > 0. When fi = . . . = /jv = x say, the latter function is 
equivalent to that previously introduced in |llj as the limiting eigenvalue density at the soft 
edge for the Chiral and Gaussian /3-Ensembles {(3 even). For general / = (/i, . . . , Jn) and 
a > 0, one has 

A W Ai H (/)=pi(/)Ai H (/)5 (5-lg) 

where A^ is the differential operator defined in Eq. (14. 28ft and p is a power sum (see Section 
2.2). This multivariate Airy differential equation follows from the use of the Calogero-like 
equation [3] 

A( Q ) ^o(x; y) = p 2 (y) Q T {x; y) (5.19) 

and simple manipulations in the integral of Eq. f)5. lTf) . 

According to Eq. (|4.22p . a multivariate Hermite polynomial associated to a rectangular 
partition is equivalent to an average of products of characteristic polynomials. By virtue of 
Proposition [TJ the latter quantity can be replaced by the expectation value of products of 
determinants in an ensemble of random matrices with an external field. Explicitly, 



xeqaEjv 

= e- s2 ((detiY)^e-^ YS > YeG/3 , En (5.20) 

if = 4/(3 = 1,2 or 4. The last integral representation suits perfectly for the large N 
asymptotic analysis. The following result generalises the well known asymptotic expansion 
of the Hermite polynomial H^{x) in terms of the Airy function. 

Proposition 10. Let Gjv,/3(S) denote the first line of Eq. (I5,20p . where S is a Hermitian 
4/ '(3 '-matrix whose eigenvalues are given by s = {s\, . . . , s n ). Then, as N — > oo ; 

^.-"'•"^(VSH-^S) . Al (^»( S ) + (_L_). (5.21) 

where 

/ A r \ nN/2 

C = 2 n(n-l)//3 ( £^ \ ^ N l/^n/2+n{n-l)/P_ ^ 22) 

The correct proof of the proposition is given in Appendix B. It is purely technical and 
relies on the multidimensional steepest descent method. However, the asymptotic formula 
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can be easily understood thanks to the following heuristic argument. First, by rescaling the 
matrix Y in Eq. (15.2Q[) . one finds 

G N ,p (V2N1 + ^=7f S ) = D f ( dY ) ex P ( N f( Y ) + N^giY, S)) (5.23) 

where 

/(Y) = tr (— 2Y 2 — 4iY + In Y) , g(Y, S) = -2tr (YS + iS) (5.24) 

and 

{ _ 2N) nN/2 {2N) n/2 + n(n-l)/p g2 

Second, the following change of variables is made in order to eliminate the quadratic terms 
in/(Y): 

W = 2iV 1/3 ^Y + il^ . (5.26) 
Third, the function /(Y) is formally expanded in powers 

f jV-i/3. This yields 

x -Afi/3 trS 1 trS 2 „ / nr= 1 \ 

J (dW) exp (I tr W 3 + i tr Ws) exp ( - g tr j . (5.27) 



Clearly, the last line is similar to the matrix Airy integral plus a correction of order less or 
equal to A r ~ 1 / 3 . However, the preceding approach is not rigorous: it doesn't proves that Eq. 
(|5.26|) furnishes the main contribution to the integral when N is large; it doesn't explain why 
W should be considered Hermitian 4//3-matrix. 

Proposition [10] only provides the dominant term to the average of products of characteris- 
tic polynomials, or equivalently to multivariate Hermite polynomial, evaluated at the edge. 
One can nevertheless find an infinite (but not convergent) asymptotic series by exploiting 
Eq. (|5,29p . Note that the latter equation remains true when considering the steepest descent 
method used in Appendix B. Recall the well known identity (see e.g. |32t Chapter I]) 

exp V -v n ) = V — v x , (5.28) 

where the sum is taken over all partitions A, v\ = v Xl ■ ■ ■ v\„, and z\ is the quantity defined 
just above Eq. (|2,10p . The use of the latter formula in Eq. (|5.29p leads to 

C L e G Nt/3 ( V2iVl + 



V2iVV3 

J (dW)exp ^-trW 3 + itrWFj ^^ yr — p A+ (W). (5.29) 

In the last equation, A + = (Ai + 3, . . . , X(, + 3) and p A (W) = tr W Al • • • trW Af . Now, let 
T stand for an arbitrary n x n matrix with no symmetry property (i.e., n 2 independent 
elements). Then 

Pk(d T )e itlWT := tr («9 T )V trWT = tr (iW)V tr WT (5.30) 
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and the next proposition follows. 

Proposition 11. With the above notation, one has formally 



c- 1 



-JVVStrS 

e 2at 



1 tr Ci 2 
Tj3 tTS Q 



f 2Nl + 



1 



E 

L A 



(-1) IAI 1 



.(^AiO^T) 



(5.31) 



T=S 



6. Conclusion 

New dualities between different ensembles of random Hermitian matrices have been ob- 
tained in the article. Their general form, valid for all (3 > 0, has been given in Propositions HJ 
ESQ and[3 

In Section 5.2, Proposition [7] has been used for proving the equivalence of the formal 
one-matrix model and the Gaussian model with an external field. But the comparison of 
Eqs. (15.201) and (15. lip leads to another surprising conclusion: the partition function (|5.9[) is 
a multivariate Her mite polynomial! Explicitly, 

(Q s t ) N Z n ,p(t u t 2 , ...)= H$X..)(°i> S2 > • • •) ^ 

i>l 

when 

/ 2\ k/2 

t k = 25 k , 2 +[±-\ pfcCsj 1 ,^ 1 ,...) (6.2) 

for all k > 1. Similarly, the partition function, for /^-matrices that are both Hermitian 
and positive (chiral ensembles), is equivalent to a multivariate Laguerre polynomial whose 
partition is rectangular. Since the partition function of a matrix model is related to the 
enumeration of maps, the above result suggests that there exists a combinatorial interpreta- 
tion for the multivariate Hermite and Laguerre polynomials. This would generalise Viennot's 
combinatorial work on classical polynomials in one variable |40j . 

Formal matrix models have been meanly studied when (3 = 2. This case is simpler and 
many tools have been developed for calculating the "large N" or topological expansion of the 
partition functions. A new and very general approach, which seems relevant to the question 
of dualities, has been introduced recently by Eynard and Orantin [19]. These authors have 
shown that the "free energy" (i.e., — InZ), for many matrix models or even algebraic curves, 
is invariant under a certain class of transformations. It would be interesting, on the one 
hand, to check if these transformations include precisely the source-external field exchange 
considered in Section 4 and, on the other hand, to determine if the invariance of the free 
energy can be generalised to the (3^2 cases. 

For j3 = 1,2,4, the large N limit of the partition function (|6.ip (equivalently, the expec- 
tation value of products of characteristic polynomials) has been evaluated in Section 5.3. It 
is proportional to a matrix Airy (or Kontsevich) integral. The proof given in Appendix B 
is valid only when one works with the matrices themselves, but not with their eigenvalues. 
However, the fact that the asymptotic formula has the same form for three distinct values 
of (3 is non trivial. Thus, it is reasonable to surmise that formula Eq. (15.2ip remains the 
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same for all (3 > if one uses Eq. (|5.17j) . for the definition of the multivariate Airy function, 
instead of Eq. (15. 14ft . The proof is still missing. 

Finally, it is remarkable that the dualities concerning the product of inverse characteris- 
tic polynomials don't involve a change in (3 (for instance, compare Propositions Q] and [2]). 
This fact certainly indicates that averages involving ratios of characteristic polynomials are 
much more difficult to calculate. Recall that the knowledge of the latter quantities gives ac- 
cess to the eigenvalue correlation functions (or marginal densities) which, from a physical or 
probabilistic point of point, are of prime importance. Note that the average of ratios of char- 
acteristic polynomials is still an open problem for the Circular /3-Ensembles |33| despite the 
fact that calculations in the latter ensembles are easier than in the Gaussian /^-Ensembles. A 
simple exercise shows, however, that the average of ratios characteristic polynomials satisfies 
a Calogero equation related to a superalgebra, gl(p\q) say, if beta is rational. Superalgebraic 
systems of the circular (or trigonometric) type have been previously studied by Sergeev an 
Veselov [411 142j . The solutions to these models are supersymmetric polynomials which belong 
to a special family of symmetric functions in two sets of variables (see examples 23 and 24 of 
Chapter 1-3 in |32j). Preliminary calculations indicate that solutions also exist for systems of 
the rational type (i.e., Calogero models connected to Gaussian /3-Ensembles) . The relation 
between the superalgebraic Calogero models and the eigenvalue correlation functions for the 
/^-Ensembles will be the subject of a forthcoming paper. 
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Appendix A. /3-Matrices 

A N x N matrix whose entries are real ((3 = 1), complex ((3 = 2), or quaternionic real 
((3 = 4:), can be represented as 

P-l 

X = X° + ^X fc e fc , (A.l) 
fc=i 

where X fc ( for < k < (3 - 1) is a N x N matrix with real elements. The "imaginary 
numbers" satisfy e\ = —1, e^ej = — e^-ej, and eie2 = 63. The conjugation is defined by 

X = X°-^X fc e fc . (A.2) 
k=l 

X is Hermitian if Xt := (X)* = X. Equivalently, X is Hermitian if 

(X )* = X°, (X. k f = -X fc k > 1. (A.3) 
The case (3 = is trivial: it is obtained when X = X° = diag(xi, . . . , xn)- 

Note that the case (3 = 4 is special since quaternions are not commutative. Thus, in general, 
Y^i j XijYj t i 7^ fo r (3 = 4. The correct definition of the trace for quaternionic 
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matrices is the following: 

i N 

tvX:=-J2(Xi,i + X i>i ) (/3 = 4). (A.4) 

i=l 

This trace is real by definition and satisfies trXY = trYX for all quaternionic matrices X 
and Y. Alternatively, one can use the representation of quaternions in terms of the Pauli 
matrices: 

trX=-trP(X) (/9 = 4) (A.5) 

where 

X°, + iXh X 2 , + iX 3 



r i N 

P(X)= P(X,i) . , P{Xij) 
L Jy=l 



~ x i,j + iX i,j x i,j ~ i x i,j 



J 



(A.6) 

In other words, -P(X) is N x N matrix with 2x2 entries. For all (3 and X Hermitian, 

detX = exp(tr lnX) (A.7) 

is equal to the product of the eigenvalues of X. Note however that, for (3 = 4, tr P(X) = 2 tr X 
and det X = v / detP(X). 

Finally, the measure on the space of Hermitian /3-matrices is simply the product of the 
real independent elements of X: 

/3-1 

!,/x! ii' /A " ini ,/A t- ( A - § ) 

i<j i<j k=l 

By using j- R dxe~ ax2+kx = ^ tt / ae k2 ^ a , one easily shows that 

m ,-.«=«XY _ (*)»»"-"" £i .Y. (A . 9) 

where both X and Y are Hermitian /3-matrices of size N. 

Appendix B. Proof of Proposition [10] 

Let S' be a fixed nxn Hermitian j3' = 4//?-matrix whose eigenvalues are close the spectrum 
edge of the G(3Ejy: 

4 = Vm + - 1 S j. (b.i) 

The aim is to determine the asymptotic behaviour of 

G^(S') = H^is') = / JJdet Ul ± Jhdj \ (B.2) 

The three cases (i.e., /3' = 1,2,4) must treated separately. Only the = 1 case will be 
detailed below since the method is easily adapted for j3' = 2 and = 4. Note that in the 
latter case, one has to use the appropriate definitions for the trace and the determinant (see 
Appendix A). 

The starting point is the rescaled matrix integral (|5.23p . The integrand is an analytic 
function in n + n(n — l)/2 real variables Y%j- By using the Cauchy theorem, integrations 
along the real axis can be transformed into contour integrals in the complex plane. By 
convention, —it < arg(Y^ < ir and the logarithm is defined on its principal branch. In order 
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to ensure the convergence of the integrals, each variable Yi j follows a path starting at ooe 1 ^ 3 
and ending at ooe 1 ^, where 

- vr < (f) s < -3vr/4 or 3vr/4 < <p s < vr and -vr/4 < (p e < vr/4 (B.3) 

The function /(Y) in Eq. (|5.23p has double saddle points, noted r/jj, if the following conditions 
are statisfied: 



° /(Y) 



= (B.4) 



r= , °' dY i3 dY M f{Y) 
together with some nonvanishing third order derivatives at the saddle point. Only constant 
solutions rjij (i.e., those that do not depend on the other rjki) are suitable for the multidimen- 
sional steepest descent method. Direct calculations imply that the saddle point conditions 
(|B.4j) are equivalent to following system of algebraic equations: 

-4r ?ii -4i + (r ? - 1 ) 4i = 0, -4 Vij + (77-%- = 0, -4 - (rj^Mv^h = (B.5) 
together with 

- 4 - (v-'Mrj-^n ~ (V'hiv^h = 0, (v^Mv^jk = (B.6) 
for all 1 < i,j,k < n and i < j. Note that these equations, 77 = is interpreted as a 

matrix. By writting the inverse in terms of cofactors, which means (r]^ 1 )^ = (det 77) cof f^j 
where det(?y) = Y^j r lij co ^Vij f° r some h ° ne shows that the unique set of solutions is: 

Vjj = ^ and rjij = (B.7) 

for all 1 < i < j < n. 

Now, one has to ensure that the descent of /(Y) is maximum when the variables Yij 
approach and leave rjij. Let 9ki denotes the argument of the complex variable Yjy = rjki and 
let 



d 3 
G1 kl 



-\ (B, 



for all 1 < k < I < n. The contours of steepest descent must comply with 

e 3ie jh +i<p jk _ _^ (B.9) 

The three possible solutions are 9^ = — 5ir/6, — ir/6, and ir/2. The two former angles are 
compatibles with conditions (1B.3|1 . Set 

W kl = 2N^ 3 (Y kl - m ) (B.10) 

Each variable Wkl follows the steepest descent path T>: it starts at ooe -57 "/ 6 , passes through 
the origin, and stops at ooe _7rl//6 . In order to get a simple expression for Gjv^S'), it is 
convenient to define the matrix W with elements W k i = Wi/-. Then, by collecting all the 
factors coming from /(Y) and g(Y,S) in Eq. (I5.23p . one shows that Gn,/3{S') is equal to 

Ce Nl/3 ^ S+ ^ s2 | p (dW)exp (-ItrW* -itrWs) cxp f - g k . k J k _ 3)/3 trW*j . (B.ll) 

where J v (dW) stands for the n + n(n — l)/2 contour integrals along the path T>, and C is 
the constant given in Eq. (|5.22p . Therefore, the following limit holds : 

lim C-V^ X/3trS_ ^^ trS2 GWS') = f (dW)exp ( --trW 3 -itrWS ) . (B.12) 

N-*00 J v \ 3 / 
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Obviously, the first correction to the above limit is of order less or equal to iV -1 / 3 . The 
complex functions involved in the integrals of the last equation are analytic, so Cauchy's 
theorem can be applied once again. The integrals remains convergent if each variable goes 
from ooe 1 ® 3 to ooe %ds where 

- vr < 6 a < -2vr/3 and -vr/3 < 9 e < 0. (B.13) 

Contours along the real line are chosen so that W can be interpreted as a real symmetric 
matrix. The change W i— > — W finishes the proof. 
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